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Abstract. We present a model for the non-thermal radio emission from bright O stars, in terms of synchrotron emission from 
wind-embedded shocks. The model is an extension of an earlier one, with an improved treatment of the cooling of relativistic 
electrons. This improvement limits the synchrotron-emitting volume to a series of fairly narrow layers behind the shocks. We 
show that the width of these layers increases with increasing wavelength, which has important consequences for the shape of 
the spectrum. We also show that the strongest shocks produce the bulk of the emission, so that the emergent radio flux can be 
adequately described as coming from a small number of shocks, or even from a single shock. 

A single shock model is completely determined by four parameters: the position of the shock, the compression ratio and velocity 
jump of the shock, and the surface magnetic field. Applying a single shock model to the 05 If star Cyg 0B2 No. 9 allows a 
good determination of the compression ratio and shock position and, to a lesser extent, the magnetic field and velocity jump. 
Our main conclusion is that strong shocks need to survive out to distances of a few hundred stellar radii. Even with multiple 
shocks, the shocks needed to explain the observed emission are stronger than predictions from time-dependent hydrodynamical 
simulations. 
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1. Introduction 



stars: mass-loss - stars: winds, outflows - radio continuum: stars - radiation mechanisms: 



About one quarter of the brightest O stars have a detectable 
radio flux of non-thermal origin, in addition to their thermal 



radio emission (Bieging et al. ll989l l. The thermal radio emis- 
sion is due to free-free emission in the stellar wind (Wright & 
Barlow [l975) . whereas the non-thermal radiation is probably 
synchrotron emission by shock-accelerated electrons (White 
[1985}. Usually, the thermal emission is only a small fraction of 
the flux, and the non-thermal emission dominates at centimetre 
wavelengths. 

In a previous paper (Van Loo et al. 120041 hereafter Paper I) 
we used a phenomenological model to explain the observed 
spectrum of the 05 If star Cyg OB2 No. 9. The synchrotron- 
emitting electrons are assumed to be accelerated by wind- 
embedded shocks, such as those produced by the instabil- 
ity of the radiative driving mechanism. A power law was as- 
sumed for the momentum distribution of the electrons and 
the synchrotron-emitting volume was assumed to extend to an 
outer boundary i?max- This approach leads to a model with 
a small number of free parameters, all of which are well- 
determined by the fit to the observations. Particularly, the in- 
terplay between non-thermal emission and thermal absorption 
places a firm constraint on the outer boundary i?max- 
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In Paper I, the spatial distribution of the relativistic elec- 
trons was given by a power law. This implies that the radiation 
is emitted by a continuous volume extending out to i?max- In 
the present paper, we take a closer look at this assumption. In 
the shock-acceleration mechanism discussed in Paper I (first 
order Fermi acceleration. Bell 1978 1 electrons are accelerated 
to relativistic energies by making many round-trips over the 
shock. Relativistic electrons are subject to a number of cool- 
ing mechanisms, most notably inverse Compton scattering in 
the presence of a dense radiation field (Chen 1992 1. They may 
therefore lose their energy rather fast when they eventually es- 
cape downstream and move away from the shock. In this case 
the emission is not expected to come from a continuous vol- 
ume, but rather from a number of (possibly thin) layers behind 
the shocks. 

Hydrodynamical models (e.g. Runacres & Owocki '2004^ 
of instability-generated shocks show a wide variety of shock 
strengths (both in velocity jump and compression ratio). We 
find that the strongest shocks generate most of the synchrotron 
emission, which means that the number of layers responsible 
for the observed synchrotron emission can actually be quite 
small. 

In the following, we first (Sect. |3 derive the momentum 
distribution of relativistic electrons in the presence of cooling 
and discuss the effect on the emergent flux. In Sect. 13 we then 
apply a model where the emission is concentrated in a smaU 
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number of thin layers to radio observations of the same star 
which was used in Paper I, Cyg 0B2 No. 9. Finally (Sect.0}, 
we discuss the implications of our results and draw some con- 
clusions. 

2. Emission layers 

2. 1. Momentum distribution at tiie stiock front 

Electrons are accelerated at a shock by the first order Fermi ac- 
celeration (Bell i 1978t . Every time an electron crosses the shock 
front, it gains a small amount of momentum. When an electron 
is kept in the vicinity of the shock front and crosses the shock 
front many times, it can be accelerated to relativistic momenta. 
The acceleration rate of the Fermi mechanism is given by 

dp Ap 

dt " 



At. 



(1) 



cycle 



where 

4 Au 

Ap-- p, (2) 

3 c 

is the momentum gain for an electron per round-trip. The 
round-trip time Atcycic can be determined as follows. After 
crossing the shock, a typical electron travels a mean free path 
A = pc/eB (White 1985) , where B is the magnetic field 
and other symbols have their usual meaning, before it is scat- 
tered by magnetic irregularities. An electron does not neces- 
sarily return to the shock after one scattering. From Lagage & 
Cesarsky ( I1983> we derive that an electron is scattered c/2u 
times before recrossing the shock, where u is the speed of the 
bulk motion upstream (ui) or downstream (U2) of the shock 
(both measured in the shock frame). Since the time between 
scattering events is given by 2A/c, the time needed for an elec- 
tron to cross the shock is X/u. Therefore the round-trip time is 
given by (Lagage & Cesarskv ll983> 



At 



cycle 



(3) 



A A 

For a shock with a compression ratio x = ui /u2 and a velocity 
jump Au — ui~ U2, the speeds (measured in the shock frame) 
are given by 

xAu 



Ul 



U2 = 



Then 
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(4) 



(5) 



eB xAit 

The round-trip time Atcycio given in Eq. is different from 
the expression Aicycie = 2A/c given by Chen (1992 ). As this 
is the expression for the mean time between scatterings, the un- 



derlying assumption in Chen J1992> is that the electron returns 
to the shock after one scattering. 

Electrons are accelerated to high momenta, but also lose 
their momentum via different energy loss mechanisms. By tak- 
ing the cooling into account, the net momentum gain rate of 



the acceleration process is drastically reduced. At the high mo- 
mentum end of the spectrum, electrons are mainly cooled by 
the inverse-Compton scattering of photospheric UV photons at 
a rate (Rybicki & Lightman [l979> 



ic 



P 



(6) 



where ctt is the Thomson cross section, the stellar luminos- 
ity and r is the distance from the star 

The net momentum gain per time interval due to first order 
Fermi acceleration and inverse-Compton cooling is then given 
by 



(7) 



IC 



The highest momentum attainable for relativistic electrons is 
determined by the balance between acceleration and inverse- 
Compton cooling at the shock. The result is a high momentum 
cut-off Pc given by 



Pc 



rricC 



xAu^ AnRl eB^ Vyot 



X - 1 CTTi* 



1/2 



(8) 



where i?* is the stellar radius, Wrot the rotational velocity of 
the star and Voo the terminal velocity of the stellar wind and 
where we used B ^ B^^^ (Weber & Davis 1967). Using 
typical hot-star values (L, = 10^ Lq, Voo = 2 500 kms^^, 
Wrot = 200 kms"\ i?, = 20 R© and Au = 100 kms^^), we 
can estimate a numerical value for the high momentum cut-off. 



15MeV / B^x 



X 



1 i?* 



(9) 



Note that p^ increases further out in the wind, due to the rapid 
depletion of UV photons. The high momentum cut-off is also 
weakly dependent on the compression ratio. Electrons can at- 
tain higher momenta in weak shocks, because the acceleration 
rate is higher in these shocks, as can be seen from Eqs. Q and 
(|5}. The high momentum cut-off is a factor of a/c/Au « 30 
lower than the value derived in Chen J 19921 ) and in Paper I, due 
to the different expressions for Atcycie (Eq. (|5|)). 

In the absence of cooling, the acceleration of electrons 
through first order Fermi acceleration results in a power law 
momentum distribution (Bell 1978 ). When cooling is taken into 
account, the momentum distribution can be quite different from 
the power law. Near the high momentum cut-off, the net mo- 
mentum gain goes to zero and the distribution will deviate from 
a pure power law. At low momenta, however, the acceleration 
rate is much higher than the cooling rate, so the momentum 
distribution is still a power law. Chen ( 1992) showed that the 
spectrum including Compton cooling can be expressed as 



N{p) = NEP-''^ ( 1 - ^ 



2\ («-3)/2 



where 

X + 2 



n = 



x-1' 



(10) 
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and Ne is a normalisation constant related to the energy den- 
sity of the relativistic electrons (see below). The shape of the 
distribution is shown by the thick solid line in Fig.^ The hook 
at the high momentum cut-off corresponds to the accumula- 
tion of electrons that, in the absence of cooling, would have 
had momenta above pc- (A downward hook corresponding to 
the depletion of very energetic electrons can also occur. For a 



detailed explanation see Van Loo (I2005> .^ As the shape of the 
momentum distribution does not depend on the round trip time, 
Eq. is unchanged from Chen (I1992t . 

We determine the normalisation constant Ne by specify- 
ing the amount of energy that goes into the relativistic elec- 



trons. Following Ellison & Reynolds (I1991I I. we assume that 
the relativistic electrons contribute little to the overall momen- 
tum balance across a shock. This means that the pressure of 
the relativistic electrons is only a small fraction {Q of the total 
pressure difference. It is unlikely that C, is more than 0.05 and it 
could be conceivably less (Blandford & Eichler fT987l Eichler 
& Usov 1993j . From the Rankine-Hugoniot relations we find 
that the pressure difference across the shock is 



AP 



Pi- 



X 



X 



1 



(12) 



Here pi is the density of the preshock gas. For simpUcity we 
assume pi — M / Airr'^Voo- The relativistic electron pressure is 
given by (White 1985) 



1. 



:C/rc 



dp pcN (p) , 



(13) 



with Urci the energy density of the relativistic electrons and 
Po the low momentum cut-off. For reasons that will become 
apparent in Sect. 13.3.31 we set C, equal to 0.05 (as opposed to 
( < 0.05). Using Eqs. il2i and J13> . we can then determine 

Ne- 

The normalisation of the momentum distribution differs 
from the one used in Paper I, which was based on the relativis- 
tic number density being smaller than the total electron number 
density. For a given momentum distribution, the total number 
of relativistic electrons and the total energy are of course di- 
rectly related. A constraint on one implies a constraint on the 
other. The energy constraint is, in principle, the more stringent 
as it is reached well before all electrons become relativistic. 

2.2. Momentum distribution beliind tine stiock front 

Once an electron leaves the shock front, the acceleration ceases 
and cooling mechanisms will reduce its momentum. This 
means that the momentum distribution of relativistic electrons 
behind the shock will differ from the distribution observed at 
the shock. To describe the distribution behind the shock, we 
need to know where the electrons were accelerated and how 
their momentum is changed under the influence of cooling. 

While the shock travels through the stellar wind with ap- 
proximately the terminal speed v^o, the relativistic electrons 
will move away from the shock front with the gas, at a signifi- 
cantly lower speed U2 (see Eq. (0}). This means that outer wind 
electrons may well have been accelerated in the inner wind. If 



the shock is at R^, the distribution of electrons at radius r, close 
to i?Si was initially accelerated when the shock was at 



|i?s 



U2 



<Rs, 



(14) 



where we have assumed that Voc and U2 are constant through- 
out the wind. This expression is valid for both forward shocks 
(which move outward through the gas) and reverse shocks 
(which move inward through the gas). The momentum distri- 
bution produced at ri is given by Eq. (I10> with pc and Ne eval- 
uated at Tj. 

Once we know where the electrons are accelerated, we 
only need to take into account the cooling of the electrons. 
Inverse-Compton cooling, adiabatic cooling, synchrotron cool- 
ing and Bremsstrahlung cooling can all play a role. Since 
the observable radio emission comes from electrons with mo- 
menta of the order ^100 MeV/c (see Paper I), we will only 
incorporate the cooling mechanisms which act at these mo- 
menta and higher. At high momenta, the electrons are mainly 
cooled by inverse-Compton cooling which dominates over syn- 
chrotron cooling (Chen 1992 1. For electrons with intermedi- 
ate momenta, the adiabatic cooling is more important than 
Bremsstrahlung (Chen [19921 . We will therefore include only 
inverse-Compton and adiabatic cooling. 

The cooUng rate of the relativistic electrons is not only de- 
pendent on the momentum of the electrons, but also changes 
with distance. Since the relativistic electrons move along with 
the thermal gas, we will express the cooling rate as the loss 
of momentum over a distance dr, where we take dr = Vocdt. 
The cooling rate for relativistic electrons is then given by (e.g. 
Chen 19921 



dp 
dr 



P 



3r' 



(15) 



where the first term represents the inverse-Compton cooling 
with flic = arL*/ {iT^rnlc^Voo) and the second term the adi- 
abatic cooling. An electron which was initially at ri with mo- 
mentum Pi will be cooled down to a momentum (Chen [T992t 



2/3 



1 



3aic 
5 



1 



5/3 



(16) 



when it reaches r. 

The above expression enables us to determine the mo- 
mentum distribution after cooling. We assume that a shock 
is spherically symmetric, i.e. that it covers a complete solid 
angle. The number of relativistic electrons initially having a 
momentum between p^ and p^ + dp^ at a radius rj is then 
AurfNijiCjdpidr, with Nijj) the differential number distribu- 
tion given by Eq. ( I10> and Aurfdr the volume that contains 
the electrons. Inverse-Compton and adiabatic losses reduce the 
momentum of the electrons to values between p and p + dp 
where p is given by Eq. (I16> . Also, the volume of the shell 
expands to Airr'^dr. Since the number of electrons in this mo- 
mentum interval and volume is not changed, the number distri- 
bution after cooling Nc [r, p) is given by 



N,{r,p)^N{pi) (^) 



2 dpi^ 

dp ' 



(17) 
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1.5 R. 

shock at 600 R. 



p (MeV/c) 

Fig. 1. The cooling of relativistic electrons behind a strong 
shock (x = 4 and Au = 100 kms^^). The thick solid line is 
the momentum distribution at the shock located at 600 i?*. The 
three solid lines are the distributions away from the shock, at 
locations 600 i?, - 1.5 i?*, 600 i?* -3 i?* and 600 i?, -4.5 i?,. 
We assumed N{1 MeV/c) = 1 at 600 and = 100 G. 
The dashed lines indicate the momentum that the relativistic 
electrons must have to produce synchrotron radiation between 
2 and 20 cm (450 MeV/c > p > 140 MeV/c). 



where pi is given by the inverse function of Eq. il6\ . Eq. (I17> is 
valid for all momenta lower than the momentum pc,e, to which 
an electron that had a momentum pc at has cooled at r. 
Therefore, pc,e is given by Eq. il6\ with pi = pc. Note that the 
momentum distribution at the shock front, N{p) in Eq. (I10> . is 
given by Nc{Rs,p)- 

An example for the momentum distribution is given in 
Fig. n] Here we plot the momentum distribution at the shock 
front for a strong shock (x = 4 and Am — 100 kms^^) 
at 600 i?* and the momentum distributions which are found 
1.5 i?,, 3.0 i?* and 4.5 i?* behind the shock. These distri- 
butions were accelerated respectively at 470 i?*, 339 i?, and 
209 i?*. It is clear that the highest attainable momentum drops 
off rather rapidly behind the shock: electrons which are 4.5 i?, 
behind the shock no longer contribute to the observable radio 
emission, i.e. the emission layer is no thicker than 4.5 i?*. 

Eqs. il6\ and (I17> are valid for both reverse and forward 
shocks. The momentum distribution at a distance \Rs — r\ be- 
hind a reverse or forward shock is initially accelerated at r^. 
Electrons with a momentum pi at are then cooled down to 
a momentum p at r. Although r is not the same for a reverse 
and forward shock, the momentum p is not significantly differ- 
ent as long as \Rs ~ r\ -C \Rs ~ Since this criterion is 
fulfilled for U2 ^ Woo, this means that a reverse and forward 
shock will have the same momentum distribution at a distance 
\Rs — r\ downstream of the shock. It is therefore possible to 
discuss only forward shocks in the rest of the paper without 
loss of generality. 



2.3. Synchrotron emission layers 

The previous section shows that the relativistic electrons (and 
thus also the synchrotron emission) are limited to narrow lay- 
ers behind the shock. The width of the synchrotron emission 
layers is determined by the distance over which the relativistic 
electrons still produce observable radio emission. The emission 
at a given frequency produced by a distribution of relativistic 
electrons is given by (Paper I) 



1 

47r 



dpNcir,p)P^{r,p), 



(18) 



where (p) is the synchrotron power of a single electron av- 
eraged over solid angle. We take into account the Razin effect 
(Razin 1960 see also Ginzburg & S vrovatskii l 1 9651 1. which is 
a plasma effect that greatly reduces the synchrotron emitting 
power. 

Although the emissivity is calculated by integrating over 
all momenta, only electrons within a narrow range of momenta 
contribute to the emissivity at a given frequency. A relativistic 
electron with momentum p emits most of its energy in a narrow 
region around the frequency (Ginzburg & S vrovatskii l 1 965> 



0.29 



3 eB / p 
An rricC \ TOcC 



(19) 



In a local magnetic field of B = 0.015 G the synchrotron emis- 
sion between 2-20 cm ( 1 5 - 1 .4 GHz) is produced by electrons 
with momenta between 140 and 450 MeV/c. The emission at 
a given wavelength is thus related to a specific momentum and 
the emission layer extends up to the distance behind the shock 
where all relativistic electrons are below this momentum. 

The width of the synchrotron emission layer ARg is deter- 
mined by the time Aicooi to cool the relativistic electrons below 
radio emitting momenta. From Eq. (I15> we see that the cooling 
rate decreases with distance from the star. This means that the 
relativistic electrons are cooled faster close to the star and that 
the width of the emission layer behind the shock is larger fur- 
ther out in the wind. The width is proportional to Atcooi. The 
cooling time A^cooi in principle also depends on the value of 
the high momentum cut-off pc, but this turns out to be a minor 
effect. Furthermore, the width of the emission layer is influ- 
enced by the outflow speed of the relativistic electrons. Because 
the relativistic electrons flow away from the shock with a ve- 
locity U2 given in Eq. 0, the width of the emission layer is 
directly proportional to U2- 

An important point is that the width of the emission lay- 
ers depends on wavelength. At short wavelengths the radia- 
tion is produced by more energetic electrons (Eq. (I19». Since 
these electrons are cooled more rapidly by inverse-Compton 
cooling (Fig.0, the emission layers at shorter wavelengths are 
narrower, as can be seen in Fig. |2] These differences change 
the shape of the synchrotron spectrum (Sect. 12. 6> . The emis- 
sion layers in Fig. |2] are somewhat larger than expected from 
Fig. [2 because the emission at a given wavelength is produced 
by electrons, not with one specific momentum, but within a 
narrow momentum range. 
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assumed po = 1 MeV/c, while the high momentum cut-off 
Pee is found by evaluating Eq. (|8j in Eq. ( I16t . 

The total emergent flux also has a contribution from free- 
free emission by thermal electrons. This contribution is cal- 
culated using the Wright & Barlow formalism (Wright & 
Barlowll975> and added to the non-thermal flux. 




r/R. 



Fig. 2. The emissivity behind a strong forward shock at r = 
600 i?, for 2 different wavelengths. The dashed line is the 
emissivity at 20 cm and the solid line at 2 cm. We normalised 
flie emissivity to the highest value. 

2.4. Emergent flux 

The emergent synchrotron flux of an electron distribution be- 
hind a shock is given by 



F 



Rs~ARs{u) 



(20) 



where D is the distance to the star, Ri, the characteristic ra- 
dius of emission (Wright & BarlowTgTs!), ARs{iy) the width 
of the emission layer at a frequency ly and G{r/R,j) a function 
that describes the influence of free-free absorption in the stellar 
wind (see Paper I). Because of the triangular shape of the syn- 
chrotron emissivity behind a shock (see Fig.|2ji, the synchrotron 
flux can be approximated, to first order, by 



(21) 



£|2 2 

All calculations are done using the exact expression Eq. j20> . 
but this approximate expression is useful in the discussions of 
the following sections. 

The synchrotron radio fluxes emitted behind a shock given 
in Eq. (I20> can be calculated numerically. In the code (which is 
similar to the one used in Paper I), the momentum integral of 
the emissivity, given in Eq. ( I18> . is computed using the trape- 
zoidal rule and the step size is refined until it reaches a desired 
fractional precision of lO^'^. The emissivity has to be calcu- 
lated at different positions behind the shock. Then, the spatial 
integral of the flux, Eq. j20> . is treated in the same way as the 
momentum integral, but up to a precision of order 1%. Here 
the function G{r/R^) is determined using a linear interpola- 
tion between precalculated values. The number distribution is 
given by Eq. ( I17> . while the mean synchrotron power of a sin- 
gle electron P^y is determined using precalculated values in the 
same way as G{r/R^). We included the Razin effect in the 
calculation of (see Paper I). The low momentum cut-off is 



2.5. Dominance of the strong shocks 

Hydrodynamical simulations show that a hot-star 
wind is filled with shocks out to very large distances 
(Runacres & Owocki l2004t . These shocks exhibit a large 
variation in shock velocity jump and compression ratio. In 
this section we will show that the shocks with high x 
Au dominate the radio emission. Consequently, the shock 
distribution responsible for the observed emission corresponds 
to a few isolated shocks in the wind. 

The normalisation constant Ne of the momentum distri- 
bution is proportional to the pressure difference AP across 
the shock, which has a Au^-dependence. This means that the 
number of electrons in all momentum ranges increases with 
the shock velocity jump. Furthermore the synchrotron emis- 
sivity jy, given by Eq. ( I18> . is directly proportional to Ne 
which means that the emissivity must be proportional to Au^. 
In Sect. l2.3l we showed that the width of the emission layer is 
also proportional to Au. From Eq. (12 1> it then follows that the 
synchrotron flux is proportional to lS.v?. This relation remains 
valid if we use the exact expression Eq. (I20> to calculate the 
synchrotron flux. 

The radio fluxes do not only increase with the shock veloc- 
ity jump, but also with the compression ratio The slope of the 
momentum distribution becomes shallower with higher x (see 
Eq. (I10» . Assuming that the number of relativistic electrons is 
the same near pq (i.e. Ne independent of x), more electrons 
are then present in the momentum range important for observ- 
able radio emission. As the emissivity at a given wavelength is 
related to the number of electrons with a specific momentum, 
more emission is radiated by the strongest shocks. A compli- 
cation arises because the normalisation constant Ne actually 
decreases with x (Sect. lzll . thereby slightly counteracting the 
effect of the slope. However, this complication turns out to have 
a minor effect. A grid of models (discussed in Sect. 13. 2> shows 
that the flux only fails to increase with increasing x if the shock 
is close to the stellar surface and if at the same time the mag- 
netic field is strong. Even in such cases the effect is limited to 
20 % on the flux. 

Thus, the shocks with the highest compression ratios and 
shock velocity jumps in the wind dominate the radio emission. 

2.6. Interpretation of the spectral index 

The broadening of the synchrotron emission layers toward 
longer wavelengths has consequences on the shape of the ob- 
served synchrotron flux: the spectral index' a of the emission 
is steeper than expected. For the sake of clarity, we neglect the 
Razin effect and the free-free absorption in the discussion. 



The spectral index is defined as Fj, oc i/" oc A 
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Table 1. Relevant stellar parameters for Cyg OB2 No. 9 
adopted in this paper. The numbers are taken from 
Bieging et al. (1989 BAC), Leitherer L) and 

Herrero et al. a999. HCVM). 



Parameter Ref. 



Viot sini (km s 
Teff (K)'' 


145 


BAC 


44 500 


HCVM 


R* (R©) 


22 


HCVM 


d (kpc) 


1.82 


BAC 


Voo (kms^-^) 


2 900 


L 


M {Mq yr-^) 


2.0 X 10"^ 


L 


L* (Lq) 


1.74 X 10^ 


HCVM 



" We follow White (1985 ( in assuming that v^ot = 250 kms~^, 
which is not in contradiction with the value for Vrot sin i. 

The wind temperature is assumed to be 0.3 times the stellar effective 
temperature (Drew fT989t . 



For a power-law momentum distribution of relativis- 
tic electrons with index n, the synchrotron emissivity 
is given by a power law in frequency with an exponent 
(Rybicki & Lightman[T979j 



Eq. i2H . with G{Rs / Rv) = 1, shows that the spectral index 
a of the synchrotron flux would equal a if the width of the 
emission layer were independent of frequency. But as the emis- 
sion layers broaden toward longer wavelengths (see Sect. 12. 3> . 
the synchrotron flux increases at these wavelengths and the ra- 
dio spectrum steepens. As a consequence the spectral index is 
lower than (1 — n)/2. 

In astrophysical applications the measured radio spec- 
tral index a is often used to estimate the compression ra- 
tio of the shock that accelerated the relativistic electrons (e.g. 
Chen & White [T994l . If we observe synchrotron emission with 
a radio spectral index a = —0.75 then for a non-cooling model 
we would derive a compression ratio x = 3 using a = a 
and Eq. il 1> . However, taking into account cooling behind the 
shock, we find that the observed flux corresponds to a strong 
shock with x = 4. 

3. Application and resuits 

3.1. Cyg 0B2 No. 9: 21 Dec 1984 data 

In Paper I we appUed the continuous model to the non-thermal 
radio emitter Cyg OB2 No. 9 (05 If). Of all the observational 
radio data that are available for Cyg OB2 No. 9, we used the 
1984 Dec 21 observation with the VLA (Bieging et al. 19891 ) 
because of its detection of emission in three radio wavelength 
bands and the small error bars on the detection. However the 
error bars on the Cyg 0B2 No. 9 observations listed by Bieging 
et al. J 19891 1 are surprisingly small (0.1 mJy). The error due to 
the flux calibration should already be 5 % of the observed flux 
for 2 cm observations and 2 % for 6 and 20 cm observations 



E 



1 1 
5 GHz 




15 GHz / 






\. 1.4 GHz 


iiiUUci pdiaiiicli;! !>■ 




X=^; Au=425 km/s; Rg=625 R. 




and B.=100 gauss 





)l 1 1 1 1 1 

5 10 15 20 25 



X (cm) 

Fig. 3. The VLA observations (21 Dec 1984) for 
Cyg OB2 No. 9 with revised error bars. The solid line 
represents a single shock model that explains the non-thermal 
radio emission. The stellar parameters are given in Table[2 

(Perley & Taylor"200Tl, which translates to 0.29, 0.15 and 0.10 
mJy at 2, 6 and 20 cm. We can therefore conclude that only the 
random errors (root-mean-square) were listed. If we include the 
calibration errors (absolute flux scale errors), we find: 5.7 ± 0.3 
mJy (2 cm), 7.4 ± 0.2 mJy (6 cm) and 4.9 ± 0.14 mJy (20 cm). 

A significant fraction of the radio emission is not ra- 
diated as synchrotron radiation by relativistic electrons, 
but is due to free-free emission from the stellar wind 
(Wright & Barlow 1975 1. Using the values from Tabled the 
contributions of the free-free emission^ to the observations of 
Cyg OB2 No. 9 are: 1.4 mJy (2 cm), 0.7 mJy (6 cm) and 0.4 
mJy (20 cm), if we assume that the stellar wind is completely 
ionised. The stellar parameters adopted are listed in Tabled We 
will now apply a layered shock model to these observations of 
Cyg OB2 No. 9. 

3.2. Single shock model 

The bulk of the emission by a set of shocks is produced by 
its strongest members (Sect. 12. 5> . To simplify the model fur- 
ther, we will consider the synchrotron emission from a single 
shock. The emission from multiple shocks can then be esti- 
mated by adding up the emission from single shock models 
(see Sect, \3.4t . A single shock model is described by four pa- 
rameters: the position of the shock Rs, the shock parameters 
Au and x which determine the momentum distribution of the 
accelerated electrons and the surface magnetic field B*. Hot 
stars are expected to have a surface magnetic field in the range 
of 10 to 200 G. The upper limit corresponds roughly to the de- 
tection limit of the field (no Zeeman splitting has been observed 
from any spectral line of any normal O star). The lower limit 

" Contrary to Paper I where we used fi — 1.4, 7 = 1 and — 1, 
we use fi = 1.3, 7 = 1.1 and — 1.3 in the Wright & Barlow 
formula. 
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600 




Fig. 4. (a) Left panel: The combinations of x, and that fit the observations for Cyg OB2 No. 9 lie in the grey area on 
the figure. Projections on two planes are plotted to situate the solutions in the parameter space, (b) Right panel: For all possible 
combinations of the model parameters that explain the radio observations of Cyg OB2 No. 9, the values for the surface magnetic 
field and the shock velocity jump Am are shown. The solid line represents lS.v?Bt^^ is constant. 



is set by the Razin effect, because otherwise no synchrotron 
radiation would be observed. 

Once a choice is made for the model parameters, we can 
calculate the synchrotron flux in the three wavelengths for a 
grid in i?s, X, Au and B^, and select those combinations of 
the parameters that fit the VLA observations within the eiTor 
bars. The emergent flux for one such model is shown in Fig.|3] 
All possible combinations of i?s, and are shown in 

Figs.|^ and|4]3. Strong constraints are produced on all four pa- 
rameters. In principle the ranges of the model parameters de- 
pend on pq, but we checked that the effect is minimal. 

3.3. Solution space 

3.3.1 . Compression ratio x and sliock position Rs 

Fig. 1^ shows that the shock parameters Rs and x are rather 
precisely determined. If the observations are to be explained 
by a single shock, then it must be strong (x — 3.5 — 4) and 
located between 520 - 640 i?*. 

The reason that we find strong shocks is as follows. The 
emission at 2 and 6 cm is barely influenced by the free-free 
absorption or the Razin effect (see below). In a simple anal- 
ysis, this means that, at these wavelengths, the exponent a of 
the synchrotron emissivity is simply the radio spectral index a. 
We can then derive the compression ratio of the shock directly 
from the radio spectrum. Using Eqs. (^J and ( I22t . we then find 
X = 1 — 3/2a. (Note that this means that weaker shocks pro- 
duce a steeper spectrum.) However, the steepening due to the 
wavelength dependence of the emission layer (see Sect l2.6> is 
not included in the above expression. This means that the com- 
pression ratio derived for a layered model will be larger than 
the compression ratio derived in the simple analysis above. 

The shock responsible for the observed emission at 20 cm 
must be located at a large distance from the star, due to the large 
free-free absorption in the wind. All radio emission emitted too 



close to the star is absorbed and cannot be observed. This ab- 
sorbing volume extends up to the radius {ly) where the radial 
free-free optical depth is unity. Using J^, Xffd?' = 1 with 
Xff the free-free absorption coefficient given by Allen J1973t . 
we find 

/ • \ 2/3 

Rr{<y) = 1.76x1028 (73,^2)'/' T-1/2 cm,(23) 

where T is in K, M in AfQ/yr, in Hz and Voa in kms~^ 
and where the other symbols have their usual meaning (Wright 
& Barlow 1975 1. Thus, the synchrotron emission must be pro- 
duced beyond i?T (20 cm) which is 450 i?*. Radiation emitted 
beyond this distance is not absorbed at 2 and 6 cm because the 
free-free opacity is smaller at shorter wavelengths. Although a 
shock is needed at large distances, it cannot be too far out in 
the wind. The free-free absorption at 20 cm would then be neg- 
ligible, resulting in a pure power-law spectrum that does not 
explain the observations in Fig.|3] 

The tight constraint on Rs presented here is only valid for 
a single shock. In Sect. 13.41 we will explain the observations 
by multiple shocks. These can be spread out over a fairly large 
geometric region, and the shock positions are then much less 
weU determined than in the case of a single shock. 

3.3.2. Surface magnetic field 

The surface magnetic field is less constrained than the com- 
pression ratio and the location of the shock. We even find fits 
for the surface magnetic fields higher than 500 G, which is well 
above the detection limit (« 200 G). The reason for this is 
the influence of the magnetic field on the synchrotron spectrum 
through the Razin effect. The Razin effect is smaU for frequen- 
cies 

1^ > 20-^, (24) 
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where rie is the number density of the thermal electrons 
(Ginzburg & S vrovatskii l 1 965> . For a strong magnetic field the 
Razin effect is negligible in the observed radio wavelengths, so 
that the synchrotron spectral shape is only determined by the 
other two parameters, i?s and x- For weaker magnetic fields, 
the Razin effect has some impact on the long wavelengths. 
Since more radiation is already taken away by the Razin ef- 
fect, less needs to be taken away by free-free absorption. This 
means that the shock must be somewhat further in the wind to 
explain the synchrotron spectrum. This can be seen from the 
bending of acceptable solutions for < 200 G in Fig.|^. 

3.3.3. Shock velocity jump Au 

Fig. shows the relation between the surface magnetic field 
i?, and the shock velocity jump Am. The shock velocity jump 
increases from 170 kms^^ (for _B, — 500 G) to 500 kms^^ 
(for = 70 G). Recall that we assumed ( = 0.05 (Sect.lO. 
For ( < 0.05, the shock velocity jump will be higher. 

The reason for the tight relation is the direct dependence 
of the emission on both Au and B^. In the absence of the 
Razin effect, the synchrotron emissivity produced by a power- 
law momentum distribution is proportional to A^sB*"^^-*^^ 
(Rybicki & Lightman I1979> . In a lower magnetic field, the 
relativistic electrons emit less synchrotron radiation and the 
number of electrons must increase (through Ne) to produce 
the same emission. Because x and Rs are nearly constant for 
all solutions (see Fig. |^), Ne can only be increased by in- 
creasing the shock velocity jump Aw. All our solutions have 
— (x + 2)/(x ~ 1) ~ 2, so to first order the synchrotron flux 
Fi, oc Au^B'^^'^ (see Fig. This largely explains the tight 
relation. The small difference for B^ < 200 G is due to the 
Razin effect. 

3.4. Multiple shocks 

The results we find for this single shock model can be re- 
interpreted in terms of multiple shocks. As in the case of a sin- 
gle shock, we assume that each shock covers a complete solid 
angle. The emission of the multiple shocks is calculated by 
adding up the emission of single shock models. From Sect. 12.51 
we know that the flux is proportional to Au^. Multiple shocks, 
each with a lower Au, can therefore also fit the fluxes provided 
the sum of their Au^ equals the single shock Au^. As a quan- 
titative example we calculate a multiple shock model with a 
surface magnetic field near the detection limit (S* = 200 G). 
The shock velocity jump required for a single x = 4 shock 
is then 290 kms^^ (see Fig. I^^). As a reasonable value for 
wind-embedded shocks is Au — 50 kms^^ (e.g. Runacres 
& Owocki r2004> . this means that we need about 200 shocks 
to produce the observed flux levels. A numerical calculation 
shows that a multiple shock model with 219 shocks (with 
X = 3.7 and Au = 50 kms^^) located between 50 i?* and 
925 i?* indeed produces the observed spectrum. 

The reason why we needed x = 3.7 for the multiple shock 
model (as opposed to x = 4 for the single shock model) is that 
free-free absorption partly counteracts the steepening effect of 



the emission layer. This is because the free-free absorption at 2 
and 6 cm is no longer negligible for emission produced close to 
the star For example, a shock near 100 i?* only contributes to 
the 2 cm flux and not to the 6 and 20 cm flux. A consequence 
is that the resulting spectrum of a multiple shock model is flat- 
ter than for a single shock model. This means that, to produce 
the observed spectrum, multiple shocks must be weaker than 
derived in the single shock model. 

The above analysis shows that we can bring down the re- 
quired velocity jump by introducing a large number (in this 
case 219) of shocks. This assumes that all shocks have equal 
strength. If the shocks have a range of strengths, the strongest 
shocks will dominate (Sect, \2.5i . and the emergent flux will 
be produced only by relatively small number of strong shocks. 
This would make it more difficult to bring the typical velocity 
jump down to a value consistent with hydrodynamical results. 
Whether the observed spectra can be reproduced by a model 
based on hydrodynamical predictions will be investigated in a 
subsequent paper (Van Loo et al. l2005> . 

For magnetic fields lower than 200 G, the shock veloc- 
ity jump derived in the single shock model is higher, hereby 
increasing the number of shocks needed to produce the ob- 
served fluxes. For B^, = 70 G the number of shocks (with 
Au — 50 kms^^) needed is about 1000. Obviously such a 
high number is irreconcilable with hydrodynamical models. 

3.5. Comparison with continuous model 

In Paper I the distribution of the relativistic electrons was 
described by a power law both in momentum and distance. 
The normalisation of the distribution was done by assuming 
that a fraction of the total number of electrons was relativis- 
tic. The solutions from Paper I that fit the observations of 
Cyg OB2 No. 9 could be divided into two groups: weak shock 
solutions with a constant relativistic electron fraction through- 
out the wind and strong shock solutions with a relativistic elec- 
tron fraction increasing outward. 

In the layered model there are no weak solutions (as can 
be seen in Fig.|4t, x > 3.5), which is an effect of the cool- 
ing on the radio spectral index. To explain this, let us consider 
the example of a multiple shock model given in Sect. 13 .41 This 
layered model fits the observations. It is also, to a good approx- 
imation, a continuous model in the sense that a large part of the 
wind is filled with relativistic electrons, with a fraction that is 
constant throughout the wind. Yet a true continuous model with 
the same parameters cannot fit the observations. The layered 
model has a steeper slope than the continuous model, because 
of the wavelength dependence of the width of the layer Thus 
the weak shock solutions that were found in the continuous 
model cannot occur in the layered model. 

For most solutions in the continuous model the fraction 
of relativistic electrons increases further out in the wind. This 
means that the bulk of the synchrotron emission is produced at 
the outside of the synchrotron emitting region as it is in a sin- 
gle shock model. The extent of the synchrotron emitting region 
Rmax in the continuous model is also comparable to the posi- 
tion of a single shock in the layered model. For 5* = 100 G, 
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^max is between 520 — 750 i?*. Therefore the continuous 
model and the layered model agree on the location of the emis- 
sion layer in the wind. 

4. Discussion and conclusions 

In this paper we refined the model from Paper I to take into 
account the cooling of the relativistic electrons and the large 
variety of shock strengths in the stellar wind. The cooling re- 
duces the emitting volume to narrow layers behind the shock. 
Since the strongest shocks dominate the synchrotron emission 
produced, only a small number of layers are responsible for the 
observed flux. We further simplified the layered model by as- 
suming that only one shock is present in the wind. Applied to 
a specific observation of Cyg OB2 No. 9, the model produces 
meaningful constraints on all parameters. The high values that 
we find for the shock velocity jump (compared to hydrodynam- 
ical models) suggest emission by multiple shocks. Note that 
there is an apparent contradiction between the large number of 
shocks needed to reduce the velocity jump and the small num- 
ber of shocks expected to produce the bulk of the emission. 

An important result of the single shock model is that the 
compression ratio (x — 3.5 — 4) and the position of the shock 
{Rs = 520 — 640 i?,) are well constrained. This means that 
strong shocks must survive up to large distances in the wind (a 
conclusion also reached by Chen 1993. This is consistent with 
the results from Paper I and from hydrodynamical calculations 
(Runacres & Owocki 2004 ). The shock velocity jump Au and 
the surface magnetic field i?, both have values which are not 
so well constrained and extend over a large range. However, 
these values are tightly related, so that the detection limit on 
i?* gives a lower limit on Au w 290 kms^^. 

The synchrotron emission is produced by the strongest 
shocks in the stellar wind. Because these shocks move through 
the wind, the radio spectrum must change in time. This may ex- 
plain the variability observed in non-thermal emitters (Bieging 
et al. In that case, the time-scales of the variability 

should be comparable to the flow time. When a single shock 
travels from 450 i?* to 550 i?*, the synchrotron emission at 
20 cm increases from to 4.4 mJy within a flow time of about 6 
days. Therefore tiie time-scale of variability in Cyg OB2 No. 9 
should be of the order of days to a week. However, it is unlikely 
that a shock would cover a complete solid angle, as it will prob- 
ably be disrupted by lateral instabilities (e.g. Rayleigh-Taylor). 
This will reduce the amplitude of variability. Furthermore, 
when the synchrotron emission is produced by multiple shocks, 
the variations due to the changing positions of the shocks could 
be completely averaged out. 

The efficiency of inverse-Compton and adiabatic cooling 
prevents electrons from carrying their energy very far from a 
shock. This however does not mean that the relativistic elec- 
trons need to be accelerated in situ. For example, electrons ac- 
celerated near 300 i?* still emit a significant amount of syn- 
chrotron radiation by the time they reach 600 i?, (see on Fig.|2] 
the emissivity at 597 i?, produced by electrons originally ac- 
celerated at 339 i?*). Thus synchrotron emission is still pro- 
duced far out in the wind even though shocks do not need to 
survive up to these distances. 



The presence of relativistic electrons in the wind also 
leads to the generation of non-thermal X-rays, due to the 
inverse-Compton scattering of photospheric UV photons (e.g. 
Chen & White 1991 1. Since there is an abundant supply of UV 
photons near the stellar surface, the non-thermal X-ray emis- 
sion is strongest in the inner wind where there is no observ- 
able radio emission. The emission shows up as a hard tail in 
the X-ray spectrum. Since the energy of the X-ray photons 
(E) is related to the energy of the incident UV photons (Eq) 
as £' « J^Eq (see Rybicki & Lightman .1979j . with 7 the 
Lorentz factor of the scattering electron, electrons with mo- 
menta between 10 — 20 MeV/c (7 « 20 — 40) are respon- 
sible for the emission between 1 — 10 keV. Since the non- 
thermal radio emission comes from electrons with momenta of 
order ^ 100 MeV/c, the non-thermal X-ray emission comes 
from less energetic electrons than the non-thermal radio emis- 
sion. Even so, the non-thermal X-ray emission is dominated by 
the strongest shocks. In this regard the conclusion by Rauw et 
al. C2002j that the compression ratio derived from the X-ray 
data is X « 1.75, is somewhat surprising. 

In order to produce non-thermal radio emission, only a 
magnetic field and relativistic electrons are required. Although 
a magnetic field and shocks to accelerate electrons should 
both be present in any early-type star, not all these stars are 
non-thermal emitters. A possible explanation is the value of 
the magnetic field. In a low magnetic field less radiation is 
produced by the relativistic electrons and the contribution of 
the synchrotron emission to the radio flux would be negligi- 
ble. This suggests that the thermal emitters have a lower sur- 
face magnetic field than the non-thermal emitters (Bieging et 
al. .l989J . 

The values found for Ait in the single shock model are 
rather high when compared to values derived from hydro- 
dynamical models, especially because these high values are 
needed far out in the stellar wind. If we re-interpret the emis- 
sion as coming from multiple shocks, the shock velocity jump 
does not have to be as high as for a single shock. For shocks in 
the outer wind, a shock velocity jump of 50 km s^^ is a reason- 
able value (e.g. Runacres & Owocki l2004l l. For = 200 G 
the observations are then produced by a reasonable number 
of shocks. In weaker fields, however, too many shocks are 
required (also, the Razin effect is stronger). Our synchrotron 
models will be confronted with hydrodynamical calculations 
in a subsequent paper. 

Both the high Au value and the location of the shock at 
600 i?* in a single shock model can be interpreted as com- 
ing from a colliding wind binary. The synchrotron emission 
is then produced by relativistic electrons accelerated, not in 
wind-embedded shocks, but in a coUiding-wind shock (Eichler 
& Usov ll993t . Even though the present model is not appro- 
priate for colliding winds, the position of the shock and the 
value of Au are close to those found in WR 147 (Dougherty 
et al. I2003t . This means that, although there is no spectro- 
scopic evidence for binarity, the observed non-thermal emis- 
sion of Cyg OB2 No. 9 is consistent with a colliding-wind bi- 
nary. For Wolf-Rayet stars, the relation between non-thermal 
radio emission and binarity is firmly established (Dougherty & 
Williams .2000i) . For O stars, the situation is less clear. About 
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half of the non-thermal O stars are apparently single. However, 
recent spectroscopic studies show that one of the archetypical 
single non-thermal O stars, Cyg OB2 No. 8A, is in fact a bi- 
nary (De Becker et al. i2004ii . Therefore, the possibility that 
Cyg OB2 No. 9 is a binary cannot be excluded. 
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